Phonon drag in ballistic quantum wires 
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The acoustic phonon-mediated drag contribution to the drag current created in the balhstic trans- 
port regime in a one-dimensional nanowire by phonons generated by a current-carrying balhstic 
channel in a nearby nanowire is calculated. The threshold of the phonon-mediated drag current 
with respect to bias or gate voltage is predicted. 
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I. INTRODUCTION 

The purpose of the present paper is to study the 
phonon contribution to the drag current in the course of 
ballistic (coUisionless) electron transport in a nanowire 
due to a ballistic driving current in an adjacent par- 
allel nanowire. Early the possibility of the Coulomb 
drag (CD) effect in the ballistic regime in quantum wires 
has been demonstrated by Gurevich, Pevzner and Fen- 
ton |l|] and has been experimentally observed by Debray 
et al. [||, g. Using the approach of Q], [g we con- 
sider two parallel ballistic quantum channels that are 
connected to two thermal reservoirs, each being in an 
independent equilibrium state. As was recently shown, 
although most of the heat from a current through the 
channel is generated in the reservoirs p] part of the heat 
is generated by the current carrying nanostructure it- 
self M via emission of phonons. Electrons penetrating 
into a biased (drive) wire from the leads are characterized 
by different chemical potentials, the situation is nonequi- 
librium and the phonons are generated by the drive wire. 

There has been done much work on the phonon contri- 
bution to the drag current in a two-dimensional electron 
gas situation S. However, we are not aware of any work 
on the phonon drag in ballistic quantum wires. The elec- 
trons in the nearby (drag) nanowire being initially in an 
equilibrium absorb the ballistic phonons emitted by the 
drive wire and the phonon drag current is created. Sim- 
ilar to the CD situation here we encounter that in the 
course of backscattering of electrons in the drive wire the 
phonons with quasimomentum hq^ = 2p„ are generated 
that in their turn are absorbed by the drag wire. The 
contribution to the current of two subbands ei + p'^/2'm 
in the drive wire and £„-|-p^/2to in the drag wire vanishes 
unless 

eV/2 > \ei-en\ 

again similar to the CD case Q. However, in the phonon 
drag we encounter the threshold 

eV/2 > sp„. 



Indeed, the energy-momentum conservation leads to 

hujq = Ep-t-hq,, — Sp 

that can be rewritten using q cos 6q = q^ as 
„ , 2ms 



\p + Tiqz I - p 



< 1 



Now taking into account that \p + hqz\,p should be in the 
vicinity of p„ ± eV/2vn we get the threshold condition 
eV/2 > spn- 



II. PHONON SPATIAL DISTRIBUTION 

We assume that the length of the nanowires L is 
much greater than the transverse dimensions of the wires. 
Therefore the spatial distribution of the emitted by the 
drive wire ballistic (nonequilibrium) phonons with wave 
vector q, SNq(r) given by the stationary Boltzmann 
equation 



sV6N„ = 7^ 



(1) 



where s = dujq/dq is the group sound velocity, tkoq is 
the phonon energy, TZ is the collision operator describing 
phonon generation, can be rewritten in the form 

TZ s 

{da: + ady)6Nq{x,y) ^ —, a= '-^ ^ tanipq. (2) 

•^x ^x 

We assume that the spatial distribution of the generated 
phonons depends only on the transverse coordinates x, y. 
The collision operator TZ can be written as 



7^ = 



-y 



V ^^ I 2TTh 

II 



'"^'V.IC.Kqx)!^ 



-^,fe(l-./^fe+r^,JiV^«] 
^&{£i,k+hq, - £i',k - fi^q), 



(3) 



where V is the volume of the channel, Cu'{q) =< 
l|e'q'"i|l' > is the matrix element for phonon induced 



transitions, Wq is the electron-phonon coupling con- 
stant that for the deformation potential interaction is 
Wq = TrA^g^/pWq, where A is the deformation potential 
constant, p is the mass density. 

According to approach in [H, p| we express the dis- 
tribution functions in Eq.(0) as the Fermi functions 
f^i^ip — fJ'^^'^'^) with shifted chemical potentials ^(^'^) = 
fizt eV/2, n is the quasi Fermi level that depends on the 
gate voltage and V is the bias voltage. 

What is concerning the collision operator spatial de- 
pendence we assume that it has nonzero values only 
within the nanowire. 

Restricting ourselves by the case T = we consider 
only the spontaneous term in Eq.(y) since there is no 
equilibrium phonons iVq' = 0. 
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2Ldfc 2 



X [/,,,.+,,^ X (1 - /,,,) - fiU^^^ (1 - rA)\ 



(4) 



Here in Eq.(y) we take into account explicitly that 
phonons are generated by the nonequilibrium electrons. 

We assume the channels have uniform cross sections, 
the origin of the system of reference being in the center 
of the current-carrying (drive) wire, the drag wire being 
shifted along the x axis by the distance D. Therefore we 
need the solution of Eq.(|) only for Sx > 0. 

The solution of the Eq. ( 2|) depends on the cross section 
geometry of the wire. Assuming, for simplicity, that the 
cross section of the wire is a circle (it is worth noting 
that the result does not change significantly for other 
geometries of the cross section) of the radius R we get 
for the coordinates outside the wire cross section 



27^ 



'^ Vi?2(l + a2)-(y-ax)2 
xe[R'^{l + a^)-{y~ax)^]. (5) 



The geometrical interpretation of this solution is physi- 
cally transparent (see Fig. |]) : let us draw a line through 
the center of the cross section of the wire having the angle 
ifs, tan If s — Oi with x axis. Now consider a line parallel 
to the already drawn one and crossing the cross section 
of the wire. The distance \AC\ from the point x^y out- 
side the cross section of the wire and lying on the second 
line to the first line is |y — axj cost/J^ = \y — ax\l\J\ -I- o?. 
The 0-function in Eq.(^ states that if this distance is 
smaller than the radius R (i.e. the second line does cross 
the cross section of the wire) the result is proportional to 
the length of the chord cut from the second line by the 
cross section, otherwise the result is zero. The similar in- 
terpretation is valid for the other geometries of the wire 
cross section. 




FIG. 1. Schematic representation of the cross section of 
the drive wire (shaded region) and the direction of phonons' 
propagation. 



III. PHONON DRAG CONTRIBUTION TO THE 
DRAG CURRENT 

Now in the drag wire due to an electron-phonon inter- 
action we have fn — fn + ^fn with A/„ satisfying the 
equation |1| 






(6) 



where v — dsn/dp is the electron velocity, /[/] is the 
electron-phonon collision term. For p > (p < 0) re- 
spectively the solution of this equation is 



A/„(z) = (z±L/2)i/[/„] 



(7) 



The boundary condition is A/[p > (p < 0)] = at 
z ~ =pL/2. The current then is given by 






^E 



dxdy 



2Ldp 



A Jo 2nh 



I[fn 



(8) 



Here the integration is over the cross section A of the 
drag wire. The electron phonon collision term is 

n' 

x{[/„'p'(l-/«p)(^p'-p + l) 
/np(,-L in'p' jjVp'— pj 

XO(£rip ~ Eti'p' + li^p'-p) 
' [J7i'p'[^ JnpJ-^'p— p' 
- fnpil ~ fn'p')(Np-p' + I)] 



xS{e. 



np ^n'p' 



hujp^p^)}. 



(9) 



In A'q and Wq we indicate explicitly only the longitudi- 
nal component of quasimomenta of phonons. The first 
and second term in Eq.(0) describes the Npi-.p phonon 
generation in the electron transition Sn'p' — > £np and ab- 
sorption of the phonon Npi-p via the electron transition 



£np —^ ^n'p' respectively. The third and forth terms de- 
scribe absorption of the phonon Np-pi (e„'p' -^ e„p) and 
generation (e„p — > Sn'p' ) of the phonon Np-pi . For the 
current induced in the drag wire by the phonons emitted 
by the drive wire we get 



n,n' 

dc[± 



dxdy f°^ 2Ldp 
Jo 



A 



2Trh 



dp' 
2Trh 



Wq |C„„/(q_L)|^ [fn'p' - fnp] 



(27r)2 

X {SNp>^pS{enp 

+6Np^p, (5(e„p - 



^n'p' I i^^p' —pj 
^n'p' '^^p—p')f • 



(10) 



Here the distribution functions /(„,«') are equihbrium 
Fermi functions f {Snp — m)- Assuming that the angle 
dependence is involved only through the phonon distri- 
bution we can take an average over the cross section of 
the drag wire and integrate over the angles ipq 



dxdy f q±dqj 



■dfqANq 



2R 



(27r)2 ^^ 'I (27r)2s7r 

pi p27T p7T/2 

X / qdqj^ TZsp / pdp / dip dipq 



(11) 
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p sm ip — ■ sm ipq 

R 



Since D/R ^ 1 we see that only small angles con- 
tribute to the integral 

— {painp - 1) < sinipq ~ pq < — (psincp + 1). 



Therefore, the result 

r.§(psinip+l) 



-%{'^~PSinip) 



dipq\ 1- psin(^ 



D 

'r 



fq 



TT R 

2D 



is proportional to the (solid) angle R/D of the cross sec- 
tion of the drag wire relative to the drive wire. This 
factor simply reflects the fact that phonons emitted by 
the drive wire should penetrate the drag wire. 



dxdy f q±dqj_ 

___ /,d,^7^.p. 



(12) 



Inserting this expression in the expression for the current 
and taking into account energy conservation laws 

OySnp ^n'p' I i^^p' —p) 
^S{eik+p'-p — Sl'k — flLOp'^p) — 

mS 



hs\p — p'\ 



jrS[k-p 



P-P 



xdi q- 



-n'p' 



-7ip 



Tis 



(13) 



that allows the integration over k we get 



J 



'8Dp^{T:shy 



E 



dp / dp' 



q'^dq 



p-p'\/h yjq^ -{P- P'/K)'^ 



\Cnn'{^.l 



X \Cw{c{l_)\^ [fn'p' -fnp] 



1 



\p-p'\ 

-^xJlp' -\-7n5 / {p—p')\^ "^ Jl'p-\-7n5 / {p—p')} 



x6[q- 



-np 



Tis 



ijlp+mf/(p—p') I, J- Jl'p'+mi/(p—p')) 



(14) 



Here S and 7 stand for (5 = £„ + e; — £„/ — £/' 7 = e„ 
£i' — £n' — £1 ■ Integrating over q we get 



J 



LmA^ 



SDp^n^ish) 



vS 



dp / dp' 



x[/n 



'p' 



In 



,Pnn'wiP,p') 



'""^ \P-P'\ 

\Jlp'-\-TnS/{p—p')\^ Jl'p-\-m5/{p—p') ) 
B {Sn'p' - £np - s\p -p'\) 
~^Jlp-\-ni''f/{p—p')\-^ Jl'p' -\-mf /{p—p') ) 
0{enp -Sn'p' - s\p-p'\)} , 



(15) 



where we introduced notation 
Pnn'ii'{p,p') = 



[£n'p' 



£ Y 

^npj 



4>n,n'iP,P') 



C/»7 



(t)n,n'iP,p'] 



Cw ( ^0n,«'(P.-P') j 



where 

4'n,n'{p,p') = yiSn'p' - SnpY - s'^{p-p') 



,^2 



(16) 



(17) 



Let us consider the case when n — n' and 1 = 1'. Then 
5 = -f = Q. 

Now taking into account that the phonons are emitted 
by the electrons having negative initial momentum [the 
nonequilibrium functions fi^p are Fermi functions with 
shifted by eV/2 > chemical potentials, i.e. fip = fP 
forp > and fip = /« forp < (/^ - f'' {e-[n-eV/2]), 
jR _ j-P(£ — [p + eV/2]))] and that the integrals vanish 
unless the Fermi functions / and 1 — / under the integral 
overlap one conclude that p' satisiying —00 < p' < —p — 
2ps, Ps = rns, contribute to the current. To avoid further 
confusion about the drag current direction note that we 



consider eV > and therefore the first term in the right 
hand side of Eq.(15) survives, otherwise if eV < only 
the second term in this equation would contribute to the 
current and due to fnp' — fnp the current would change 
the sign. We have after transformation p' -^ —p' 



J 



8Dp^TT^{sh) 



tE 



dp 



dp' 



^[jnp' Jnp\ I _ /I //p'(l hp) 



LmA^ 



8Dp^T:^{sh) 



vT. 



nl 



dp I dp' 

Jp+2ps 

x[Q{p.^^p')-Q{p^-p)] 
Q{p-Pi )Q{p{ -p), 



(18) 



p + p 
where 

Pn = \/2™(/i-e„), pt == V2m(^ - £„ ± eV/2). 
We get the nonzero result only if 

p+ - p- > 2ps, Pi^ < Pn < pi' (19) 

The last inequality is equivalent to 

eV/2 > |e„-e;| = |£„;|. 

Assuming eV/2 <^ p'^/2m, |e„;| <C p^/2m we have 
nonzero result only if 

eV/2 > max{sp„, |e„i|}. 

Let us put Eni = 0, we obtain for eV/2 < 2sp„ {a = 
eV/ispn), i.e. 1/2 < a < 1 



J = Jo5](— )'/ dt{2a^t)T{t) 
and for eV/2 > 2sp„ {a > 1) 



(20) 



+ / dt{2a - t) \ T{t) 



Jo 



T{t)^ 



f" 



VW^T 






2pn 



(21) 



(22) 



Cnn [^V^^ 



(23) 




FIG. 2. Drive voltage dependence of the phonon contribu- 
tion to the drag current for different values of p„. 



Assuming the following model dependence for C„ 

1 



\Cnn{qr = 



1 + qRy 



(24) 



we plot in Fig. the drag current versus the voltage 
applied across the drive nanowire for different p„ val- 
ues piR/h — 3.33 and p2R/fi = 5. Near the threshold 
a — 1/2 ^ 1 we get, assuming 2a — 1 <C (hl2pnRY 
so that the argument of the function C„„ is small and 

C n^ 1 



J = 



4x/2Lto5A 



iDpH^h 



H?) 



/ eV 



\2spr. 

1 < 



- 1 

eV 

2spn 



3/2 



(<2) 



(25) 



i.e. at the threshold the current increases nonlinear ly 
with the applied voltage. This dependence is illustrated 
by the small inlet in Fig. since it can be noted only in 
very small vicinity near the threshold. 

The second derivative of the current with respect to 
the bias voltage diverges as 



dV^ 



1 



^/eV/2spn 



1 



at the threshold for any new subband n, this fact can be 
instrumental for the experimental investigation. 

Assuming the following parameters eV/4:spn > 1, (for 
s ~ 3- 10^ cm/s, U/pn — 0.5- 10~^ cm this means volt- 
ages V ~ mV), PuR/Tt- ~ 1, L ~ 2 ^m, D ^ 0.1 • ^m, 
A ^ 8 eV, m = O.OTtoq we get the following estima- 
tion for the contribution of any subband to the phonon- 
mediated drag current 

J- 3- 10"^^ A 



IV. CONCLUSION 

We note two essential differences between the phonon 
drag and the Coulomb drag. First, it is the existence of 
the threshold eV/2 > spn (that can be achieved chang- 
ing either the bias voltage or the gate voltage). Second, 



the weak dependence on the distance D between the cen- 
ters of the drive and drag wires rather than the exponen- 
tial one in the CD case. On the other hand, note, that 
the distance dependence of the phonon drag current in 
our case is stronger than the distance dependence of the 
phonon drag between two 2DEG layers ^ . 

The coupling coefficient for the piezoelectric coupling 
for GaAs having the cubic Td symmetry can be written 
as 



W^^ 



pWq 



47re/?' 



{F{d,,V,)f 



where F{9q, (pg) is a function of angles. The estimations 
for GaAs (/3 = 10^ SGS units, e = 12, A = 8 eV, p = 5 
g/cm'^, s = 3 X 10^ cm/s) shows that the piezoelectric \°i 
interaction is the dominating one for frequencies Wq < 
10" s~\ 

Frequencies of transmitted phonons in the phonon drag 
are ojq > 2s/{h/pn), i.e. are greater than 10^^ s~^. 
Therefore, the contribution of the piezoelectric coupling 
is smaller than the considered deformation potential cou- 
pling. 

And finally, one can express the drag current also as 
an integral over the transferred phonon momentum q^ 
in an electron-electron interaction. We get a physically 
transparent result 

1 f dxdy 



J = e- 



V2 



A 



E 



w„ 



n,/,q,g3>0 



X |Cn(qi 



|2 A 
sqj_ 



i?2 



'{yqx/q-xqy/qY 



X |C„„(q±)|^e(ey- hLL>q)xb{£n, £n,lT-(^q, qz) 

xxa(ei - eV/2, ei + eV/2, fiuj^ - eV, q,) 



where 



Xai£a,£t3,'huj,q) = 



is 



E 



/(£ap) 



huj + (Eqp) - (Sap-Hq) - iO 



(26) 



(27) 



0<p<hq 

is the Lindhard susceptibility subject to the constraint 
< p < hq, while in the Xb there is no the second con- 
straint and the sum is over p > 0. This form clearly 
demonstrates that the drag current is a convolution of 
the spontaneous polarizations within each quantum wire. 
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